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J.D.  Clayton 

In  the  context  of  finite  deformation  mechanics,  a  tangent  mapping  is  anholonomic 
over  some  domain  when  it  is  not  a  gradient  of  a  motion;  conversely,  a  deformation 
gradient  is  holonomic  when  it  is  integrable  to  a  motion  field  everywhere  in  that 
domain.  This  brief  report  addresses  covariant  differentiation  for  four  possible 
choices  of  basis  vectors  in  anholonomic  space.  As  an  example  from  continuum 
physics,  the  theory  is  applied  towards  description  of  divergence  of  the  heat  flux. 
An  extensive  treatment  of  anholonomic  mathematics  can  be  found  in  a  recent 
article  [1];  however,  this  report  includes  material  not  found  in  [1],  and  vice-versa. 

As  suggested  by  Schouten  [2],  consideration  of  differential  geometry  of  anholo¬ 
nomic  spaces  dates  back  to  at  least  1926  [3].  Many  important  identities  are  derived 
in  [2,  4],  Various  coordinate  systems  and  associated  metric  tensors  in  anholonomic 
space  are  considered  in  [5],  with  particular  focus  on  convected  basis  and  Cartesian 
representations.  Correspondences  among  mathematical  objects  from  differential 
geometry  and  their  continuum  physical  counterparts  in  defect  field  theory  of  crys¬ 
tals  are  described  at  length  in  a  more  recent  monograph  [6]. 

The  present  description  is  limited  to  the  time-independent  case,  such  that  spa¬ 
tial  coordinates  xa  are  related  to  reference  coordinates  XA  by  one-to-one  and  at 
least  twice-differentiable  mappings  xa(X)  and  XA(x),  with  X  a  material  particle 
and  x  its  spatial  representation.  Let  the  usual  holonomic  deformation  gradient 
F(X)  be  decomposed  multiplicatively  as 

(1)  F  =  F  F ,  F\  =  iM; 

(2)  F  =  dAxaga  (g>  Ga,  F  =  F?aga®~ga,  F  =  F?A~ga®GA. 
Denoting  dA  =  d/dXA  and  da  =  d/dxa,  partial  differentiation  proceeds  as 

(3)  aa(-)d=  F?M)  =  F~1AdA(-),  dA(-)  =  dAxada{-)  =  F?M). 

Attention  is  restricted  to  a  simply  connected  domain  in  reference  and  current 
configurations  such  that  {XA}  and  {xa}  are  global  coordinate  charts.  Let  indices 
in  brackets  be  skew,  e.g.,  A[AB]  =f  \{Aab  —  ABa)-  Since  XA  and  xa  are  holonomic 
coordinates, 

(4)  d[AdB](-)  =  0,  d[adb](-)  =  0-,  d[AF»B]  =  0,  d[aF-\A=  0. 

Similar  identities  do  not  always  hold  for  da(-)  since  F  1  and  F  are  not  necessarily 
integrable  functions  of  xa  or  XA.  Anholonomic  object  k  obeys  [1,  2] 

(5)  K&  =  -F-\adwKx]  =  -F^F-\A  = 

(6)  d[ad, 3](-)  = 

Holonomic  charts  { xa(X )}  [or  {ia(a:)}]  exist  in  a  one-to-one  fashion  with  X  ox  x 
if  and  only  if  =  0  throughout  a  simply  connected  domain. 


Let  A  be  a  generic  differentiable  tensor  field.  Covariant  differentiation  in  an- 
holonomic  space  is  defined  as 


V,,A° 


■  ■<t>  d£.f  O  Aa. 

..fi  ~ 


i  AP- 

.p  '  X  Up  ^7- 


..p'  '  up^^.-.p  x  u-y^p.. 


(7)  ._7 

Connection  coefficients  can  be  expressed  in  general  form  as  [2] 

(8)  =  ^a,5(9{/3^x}  -  2T{0Sx}  +  2k{/35x}  +  M{0Sx}), 

where  gaxgX(3  =  Sp  and  the  following  definitions  apply: 
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(10)  9a(3  —  9cx  '  9(3)  ^bx  ~  ^[Px]  K(^X')  ^aPx  ~  ~^a9pX- 

In  this  report  attention  is  restricted  to  metric  connections  so  that  Mapx  =  0  and 

(11)  d*90  =  K}9x,  dagfj  = -r^gx- 

(12)  da  In  sjl  =  Tmf0,  VaepxS  =  dae0xS  -  F^e0xS  =  0; 

(13)  g  d=  det (ga0),  ea0x  d=  \Tgea0x,  ea0x  d=  {l/^)ea0x. 

The  covariant  derivative  of  a  generic  differentiable  vector  field  V  =  Vaga  is  then 

(14)  V  V  =  dp  V  ®  g0  =  (d0Va  +  r 0axVx)ga  ®  ~g0 . 

Total  covariant  derivatives  of  two-point  tangent  mappings  F  and  F  1  are  [1,  6] 
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(16)  F-\«a:b  = ehF-'Z  +  TftF-  mb,  .a-,ba,  mC-f  ma:p.  mb. 

Metrics  and  Levi-Civita  connections  in  reference  and  current  configurations  are 


(17) 


dcf 


Gab  =  Ga-Gb  =  dAX  ■  dBX, 


p..A  dcf  l  /-iADf)  /~i 
1  BC  ~  °{B^rDC}’> 


(18)  gab  =  ga  -  9b  =  9ax  ■  dbx ,  rb°  =  \gadd{bgdc}. 

Letting  g  =  det (gab)  and  G  =  det(GJ4s),  Jacobian  determinants  are  [5,  6] 

(19)  J=^g/Gdet(dAxa)  =  Jj,  J  =  ^gJ^det(F^),  J  =  J§/Gdet{F%). 
Piola’s  identities  for  possibly  anholonomic  Jacobian  determinants  are  then  [1,  4,  6] 

(20)  (JF-'iU  =  ea0xeABCF0AFxB:C],  (J"1^  ):o  =  ea0xeabcF-10F~lx:c]. 

Let  q  =  qaga  denote  the  heat  flux  vector  referred  to  anholonomic  space,  let  ka0 
denote  a  covariant  constant  positive  semi-definite  tensor  of  thermal  conductivity 
with  the  particular  form  dictated  by  the  material  symmetry  group,  and  let  9  denote 
temperature.  Nanson’s  formula  and  energy  invariance  among  configurations  lead 
to  relationships  among  q,  spatial  flux  q,  and  reference  flux  Q: 

(21)  qa  =  jF~\Zqa  =  J~1FaQa  =  -ka0d09. 


899 


Heat  transfer  per  unit  anholonomic  volume  is  the  divergence  [6,  7] 

(22)  +  qa' + 

=  J~1\7aQa  =  J\7aqa. 

Four  choices  of  basis  {ga}  are  considered.  In  the  first  case,  the  anholonomic 
object  is  assumed  to  vanish  such  that  Euclidean  position  vector  x(X)  exists: 

(23)  gQ  =  dax,  r  0*=l~gaSd{/3~gSx},  {JF~1±)-A=  0,  =  0; 

(24)  V«r  =  8aqa  +  qada  In  ^  =  -ka/3(dad06  -  T^dx9). 

In  the  second  case,  Cartesian  intermediate  bases  {ea}  are  assigned,  but  tangent 
maps  need  not  be  integrable: 

(25)  ga  =  ea ,  ga0  =  Sa0,  =  0,  Va(-)  =  8a(-); 

(26)  Vaqa  =  8aqa  +  qaJda(J~1F-a)  =  -k^[dad09  +  Jda{  J^F^dpO]. 

In  the  third  case,  {ga}  are  chosen  coincident  with  reference  basis  vectors  {Ga}', 
object  torsion  Xjy  “ ,  and  curvature  from  Pg“  all  may  be  nonzero  [1];  and 

(27)  ga  =f  6^ Ga,  T^  =  F~\p^r^,  VaV*  =  F~\iV  AVB  S*-, 

(28)  \7aqa  =  -ka0[dad0e  -  F-\i6%6*TACBdxe  +  j-\jF-x±):Ad0e}. 

In  the  fourth  case,  {ffQ}  are  chosen  coincident  with  spatial  basis  vectors  {ga}\ 
object  torsion  Xjy  “ ,  and  curvature  from  all  may  be  nonzero  [1];  and 

(29)  ga  =  6aaga,  r  =  VavP  =  F?aVaVbsP-, 

(30)  \7aqa  =  -k^[dad0e  -  K6b6^bdxe  +  JiJ^F^dgO}. 

The  second  case  (Cartesian)  is  most  common  and  presumably  most  practical  for 
materials  of  arbitrary  anisotropy;  the  third  or  fourth  cases  may  prove  useful  for 
structures  with  curved  shapes  and  hexagonal  or  isotropic  symmetry. 
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